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MATTHEW CARGO 

Abstract. We introduce a method for calculating the joint spectra of func- 
tions which comprise a quantum integrable system under a deformation quan- 
tization star product. The main result involves a construction by formal power 
series in h of s'Darboux coordinates, a concept we introduce for the star prod- 
uct analogue of Darboux normal coordinates. We will find that underlying 
the quantum integrable system is a set of s'Darboux coordinates which in turn 
give rise to number operators for the system. We present an explicit correction 
to the lowest order Bohr-Sommerfeld or EBK quantization rule. 



1. Introduction 

Let us describe a path from canonical quantization to deformation quantization. 
T*(R) is a symplectic manifold with natural coordinates (x,p) and Poisson bracket 
is {x,p} = 1. To quantize, we promote (x,p) to operators (x,p) on L 2 (R) with the 
canonical quantization relation [x,p] = ih. One such choice has p = —ihd/dx. We 
may now build operators out of the generators (x,p). Classical mechanics is left 
behind — only the names of operators (x,p) betray their classical origin. 

However, it is possible to recover the original classical structure by using a symbol 
correspondence, which associates to each operator A its symbol A, a function on 
classical phase space taking values in power series in h. The operator product AB 
becomes a non-commutative product of functions A-kB, called a star product. When 
A and B are suitably well-behaved, a star product may be developed as a power 
series in H of bi-differential operators. All symbol correspondences are such that the 
associated star commutator has {A, B}+ = A* B — B * A = ih{A, B} + 0(h 2 ) — the 
Poisson bracket has reappeared. 

Having thus returned to the classical phase space, it is natural to ask whether 
the Hilbert space L 2 (R) was needed at all. So began the so-called deformation 
quantization program of 1 . On an arbitrary symplectic manifold, not necessary 
a cotangent bundle, they sought to create a star product from scratch, assuming 
only that 

• the star product is associative; 

• the star product is the ordinary product at lowest order, i.e. A * B = 
A-B + 0(h); 

• the star commutator is the Poisson bracket at lowest order, i.e. {A, B}+ = 
ih{A, B} + 0(h 3 ) 

Since the inception of deformation quantization, investigators have affirmed the 
existence of a star product. We are particularly attracted to the construction of 
Fedosov 0, who builds a star product from the symplectic form and a symplectic 
connection. 



2 



MATTHEW CARGO 



In view of these successes, one might be tempted to conclude that the deforma- 
tion quantization program is complete. We disagree. As a practical matter, many 
calculations that are straightforward in traditional Hilbert space quantum mechan- 
ics are not yet possible in star product quantum mechanics. In particular, we lack 
systematic methods for computing spectra. 1 

This paper is a step in that direction. In certain types of quantum integrable 
systems, we will show that the joint spectra can be developed as a formal power 
series in h. 

The organization is as follows. In Section [31 we will introduce s'Darboux coor- 
dinates, the star product analogue of Darboux normal coordinates, and show how 
s'Darboux coordinates may be developed from Darboux coordinates. In Sections 01 
and |3 we will review Dirac's creation and annihilation operator method and define 
quantum integrable systems. In Sectional we will show that creation and annihila- 
tion operators underlie some quantum integrable systems and develop methods for 
determining these operators. In Section [7| we will demonstrate these techniques 
with the Moyal star product and show that they lead to higher order EBK quanti- 
zation rules. Finally, in Section|Hl we will generalize this result to the Fedosov star 
product. 

2. Definitions and conventions 

Symbol correspondence. Operators will wear hats, e.g. A, and we denote the 
corresponding symbol by A — s(A). When the symbol A has a power series of the 
form A = a + hA\ + ... we call the first term, or principal symbol, it A. When we 
need to name a principal symbol, we will try to use the corresponding lower case 
letter, as above where a = it A. 

Lists. On our 2M dimensional symplectic manifold, we shall often work with 
lists of M or 2M functions. When we refer to a specific item in a list, we shall 
attach an index, e.g. H l , but when we refer to the items collectively, we will drop 
the index, e.g. H . 

Other notations will be developed as needed. 

3. From Darboux to s'Darboux coordinates 

On a 2M dimensional symplectic manifold M. , Darboux normal coordinates are 
coordinates z a satisfying 

(3.1) {z a ,z h } = J a \ 
where J ab is the constant matrix 

(3.2) J ab =( T ° ImxM \ 

\ —-IMXM U J 

A collection of 2M symbols Z % are said to be s'Darboux coordinates if they have 



(3.3) {Z a ,Z b }+ = ihJ ab , 

^The present form of this paper is a somewhat rushed presentation of results which should 
have been published much earlier. I present them now to establish whatever priority they may 
have. Based on a preliminary survey of the literature, I believe the results are new. However, I 
am aware that many others have also pursued work in this direction. In a future version intended 
for official publication, I fully intend to give due credit to earlier work — hoping, of course, that I 
do not find the exact result elsewhere! 
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that is, modulo the ih, they satisfy the Darboux coordinate algebra under the star 
bracket. 

In this section we prove 

Theorem 1. For every set of 2M Darboux coordinates z l £ C°°(Ai) there exist 
s'Darboux coordinates Z l £ C°°(.M)[[/i]] such that z 1 = TtZ 1 . The Z l are not unique: 
instead, given s'Darboux coordinates Z l there are other sets Z H = Z l + h n X l with 
ttX 1 = {7, z 1 } for some 7 £ C°°(M). 

The proof is by induction; we construct Z l order by order in H. Already, 

(3.4) {z\zi}* = ihJ ij + 0(h 2 ) 
since z l are Darboux coordinates. We next assume that 

(3.5) {Z\ Zi}* = ihJ ij + 0(h n+2 ). 

To complete the induction, we need to show there exist X 1 £ C°° (Ai)[[h]\ so that 

(3.6) {z i + h n+1 x\ z j + h n+1 x j }* = ihj iJ + o(h n+3 ). 

(We then take Z l Z i + K n+1 X\) Rewriting the left side of Eq. (E3J), we have 
(3.7) 

{z l + h n+1 x\ z j + h n+1 x j }* = {z\ z j }* + h n+1 {{x\ z j }* + {z\xi}*) + o{h 2n+3 ) 

= {Z\ Z j }* + ih n+2 ({nX\ Z J } + {Z\irX j }) + 0(h n+3 ) 

If Eq. (|3.6|l is to hold, the terms containing irX 1 must cancel off the 0(h n+2 ) part 
of {Z\P}*. 

The condition can be expressed more geometrically. First define coordinates 
Zi — JijZ° using the inverse J+j of J v . Then define a two form, 

(3.8) £o n = ^-{Z\ Z 1 }* dzi A dzj 
and a one form 

(3.9) 9 = X l dzi. 
Then, since 



(3.10) d9 = P k y 1 dzj A dzi = {z\ X 1 } d Zj A dz { 



.dX, 

Eq. H3.6f) requires that 

fcra+1 

(3.11) Q + -— d6 = u + 0{h n+2 ), 

where lo = J lJ ' dz{ A dzj is the usual symplectic two form. Taking the exterior 
derivative, we obtain the condition 

(3.12) du> = ^-{z k , {Z\ Z j }*}dz k A dz, A dzj = 0(H n+2 ). 

Now, if the right side of this expression had not a Poisson bracket but a * bracket, 
and if also z k were Z k , then it would vanish by the Jacobi identity. To investigate 
whether these changes leave it unchanged at this order, we form the difference 

(3.13) 

{z k , {Z\ Zi}*} - ±-{Z k , {Z\ Zi}*}* = {z k , {Z\ Zi}*} - {Z k , {Z\ Zi}*} + 0(h n+3 ) 
in 

= h n+3 , 
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where we have used the induction hypothesis of Eq. I|H.5|I . Using this calculation, 
we have that 

diu = — {z k , {Z\ Z J }A dz k A dz t A dzj 
in 

( 3 - 14 ) = 77^2 {Z k , {Z\ Z>}*}*dz k A dz t A dzj + 0(h n+2 ) 

= 0(h n+2 ), 

where, as planned, we have used the Jacobi identity. Now, because we have assumed 
that the Darboux coordinates z % are C°° , we have also assumed that the cohomology 
of M. is trivial. Thus, u> is not only closed but exact at 0(h n+2 ), and there exists 
a 9 making Eq. l|3.11fl true. 2 Of course, 9 is defined only up to an exact one-form; 
adding d 7 to 9 adds h n+1 {z\j} to Z\ 

4. Digression: Dirac algebra 

Suppose we have a set of M operators A % which act on the Hilbert space L 2 (R M ) 
and satisfy the algebra 



(4.1) 



A l \ti 
A\A j 



= hdij 



0. 



The M operators N l = A^A 1 are positive definite and commute. Dirac showed 
the N l have simultaneous eigenspaces V n labeled by non-negative integers n % : for 
4> € V n , N l ip = hn l ij). 

We can attempt to duplicate this construction using our s 'Darboux symbols Z l . 
We create 

(4.2) A i — -^-=(Z i — iZ i+M ) 

V2 

for i = 1, . . . , M. These symbols satisfy Dirac's algebra under the star bracket: 

(4.3) 1 . / 3 

{A\A>}* = Q. 

We can also define M commuting symbols N l — A 1 * A 1 . 3 

In order to conclude that the spectra of the iV l are bounded from below, Dirac's 
original argument makes use of the positive definitcness of the N l . In this case, 
because we do not have a Hilbert space, we cannot complete Dirac's argument. 
Instead, we will simply postulate the following: if a set of M symbols A 1 £ 
C°°(M)[[h}} satisfy Eq. (|4.3|) . then the spectrum of h~ 1 N l consists of all non- 
negative integers. 4 



One could also work on a chart on Ai. We would then have a theorem about local s'Darboux 
coordinates on a manifold with arbitrary cohomology. 

3 There is no sum here. However, given a basis a a for the Lie algebra u(M), we may define 
o~1jA x * A 3 , a representation of u(M) under the * product. The N l are a maximally commuting 
subalgebra. 

^Readers are encouraged to explore if this postulate is necessary. 



S'DARBOUX COORDINATES AND WKB APPROXIMATIONS IN DEFORMATION QUANTIZATION 



5. Quantum integrable systems 

In classical mechanics on a 2M dimensional symplectic manifold, an integrable 
system is a collection of M Poisson commuting, functionally independent functions 
h l . We will define a quantum integrable system as a collection of M star commuting, 
functionally independent functions H % € C°°(.M)[[fi]]. 

Because {£P,iP}± = implies {irH 1 ,nH 3 } = 0, we may associate to a quan- 
tum integrable system H l the classical integrable system h l — irH 1 . Treatment of 
the quantum problem may begin only after understanding the underlying classical 
mechanics. 

In this paper, we will concentrate on the case where the classical evolution takes 
place on tori. To be precise, we have the usual energy-momentum map h : M. — > 
M. M : p i— > h{p). For values E in the interior of hA4, we will assume that hr x E is 
always a single M -torus. That is, the h % give rise to a global set of action-angle 
variables (0\ . . . , 9 M , I 1 , . . . , I M ) satisfying 5 

{9 l ,9 3 } = 

(5.1) {P,P} = 

{9\P} = S ij , 

and functions p : R M — » K which express the h % in terms of the P: 

(5.2) h l = P o I. 

In view of Eq. I|5.1|) . it is tempting to regard the (6,1) as Darboux coordinates. 
This we should not do, however, as the 9 l are not C°°. Nevertheless, we may create 
Darboux coordinates from this set by defining 



(5.3) 



it is straightforward to verify {z l , z 3 } — J lJ . In doing so we must assume that the 
basis contours used to define the P have been chosen so that P > everywhere. 
In other work, my collaborators will argue that such contours do exist for this type 
of integrable system. In addition, we shall assume the 9 l can be and indeed arc 
chosen so that the z % are C°°. 6 

Previously, we proved Darboux coordinates z l can be quantized to s'Darboux 
symbols Z % satisfying z l = -kZ % . In Eq. I|4.2|l and following, we showed how to 
create number operators N z from the Z % . Next, we will show that the Z % can be 
made so that {fP, JV J '}* = 0; when this condition holds, we shall say the N % are 
good number operators for H l and that the Z % are compatible with H l . 7 

6. Existence of good number operators 

In this section, we prove two theorems. The first concerns "approximately good" 
number operators. It will help prove the second, our main result. 

^In we have an exception to using the lower case for classical functions. 

^The angle coordinates l are, of course, not unique. Taking 6 1 — > (9 1 + {0 1 , 7} for some function 
7 gives new angle coordinates. 

^Already, it is possible to construct quantum integrable systems from a classical one: from h 1 
find (1,0), form z 1 , quantize z l to Z 1 , form number symbols N l , and then define H 1 = h' o N. 
Because the Z % are not unique, neither are the resulting H r . 
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We first introduce some new notation. Let / be an analytic function with Taylor 
series f(x) = ^c n x n . Define / o A = ^c„A*™; that is, for each term x n , we 
substitute the n-fold star product of A. If A is the image of A under a symbol 
correspondence, then / o A = s(f(A)). Thus the o notation gives us a way to 
talk about the symbol of a function of an operator in the context of deformation 
quantization, where there is neither a Hilbert space nor operators on it. 

For the purposes of this section, the only fact we need to know about o is that for 
★ commuting functions N l , FoN = FoN + 0(h). This fact follows immediately 
from the deformation quantization postulates. 

The two theorems can now be stated as 

Theorem 2. Suppose H l and N* satisfy {H\ W}„ = 0, {N\ W}* = 0, {N\W} ± = 
0(h n+2 ), and there are f l such that ttH 1 = f l o nN. Then there exist F l £ 
C°°(R M )P] andG i G C°°(M)[[h}} such that H* = F l o N+h n+1 G l . Furthermore, 
nG i = {G, h 1 } + g i o I, for some G,g l € C°°(M). 

Theorem 3. There exist s'Darboux symbols Z % compatible with the quantum inte- 
grable system H l . 

Before proving these we remark that nothing in Theorem [2] requires that the N l 
are number symbols. However, if they are, the spectrum of H l — F l o N + 0(h n ) 
is just F l o (hn) + 0(h n ), where are non-negative integers. 

6.1. Proof of Theorem EJ We prove the theorem by induction, constructing F z 
order by order. We begin with F l — f % . Then 

(6.1) F l o N = f o N + 0(h) = f o irN + 0(h) = H l + 0(h). 
We next assume that F l is such that 

(6.2) F l o N = H l + 0(h s+1 ). 

To complete the induction, we need to show there exist X 1 E C°°(M. M ) such that 

(6.3) (F i + h^x 1 ) Zn = h 1 + 0(h s+2 ) 

whenever s < n— 1. (We then take F l — > F l + h s+1 X % .) Working with the left side 
of Eq. Hfi.^fl , the condition becomes 

(6.4) F l oN + h^X 1 o ttN = W + 0(h s+2 ). 

Because X l onN depends only on the TV 1 , X 1 may be used to enforce this equation 
only when the other terms also depend only on the N*. To verify this, we compute 

{-KN\F i oN-H i } = -^{N j ,F l oN- H% + 0(h s+2 ) 
ih 

( 6 - 5 ) = -i (o + 0(h n + 2 )) + o(h s+2 ) 

ih 

= o(h n+1 ) + o(h s+2 ). 

This is of 0(h s+2 ) whenever s < n — 1. At the final stage of the induction process, 
we have s = n — 1 and 

(6.6) iT = F*o N + h n+l G\ 

for some G\ 
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This concludes the first part of theorem [3 For the second part, we compute 
Q = {H\H 3 }* 
= {F l *oN + K n+1 G\F 3 oN + h n+1 G j }± 

(6.7) 

= h n+1 ({F l o N, G ] }± + {G\ F 3 o TV}*) + 0(h 2n+2 ) 
= ih n+2 ({f o itN,itG 3 } + {ttG\ f o TriV}) + 0(h n+3 ). 

Thus, since {<kH\ TrG j } + {ttG\ttH 3 } = 0, we have nG 3 = {G,ttH 3 } + gi o nl 
for some functions G and g 3 . Note that G is not determined uniquely in this 
decomposition — we could add a function of the irN 1 . With this result, equation 
(|6.6I) becomes 

(6.8) H l = F l o N + h n+1 {{G, irH 1 } + g 3 o ttN) + 0(h n+2 ). 
This completes theorem 13 

6.2. Proof of Theorem |3J We will construct the symbols Z l order by order in H, 
taking care that the associated number symbols N l commute with the H l at the 
current order. We assume that Z % are s'Darboux symbols, but that {A 1 * A 1 , if 3 }* = 
0(h n+2 ). For the induction step, we will create new d'Darboux symbols 

(6.9) Z H = Z i + h n+1 g\ 

so that the associated N H = A' 1 * A H satisfy {A H * A'\ H 3 }* = 0{h n+3 ). 

As discussed before, the the freedom in the Z l , allows us to take ng l = {7, z 1 }, 
with 7 a real function. This modifies the number operators at order K n+ . First, 

(6.10) A H = (Z H - iZ H+M ) = A i + ft" +1 {7, ttA 1 } + 0{h n+2 ). 

v2 

Then 

N H = {A 1 + K n+1 {j, irA 1 } + 0{K n+2 )) * (A 1 + h n+1 {j, irA 1 } + 0{h n+2 )) 

(6.11) = N l + h n+1 (^{7, ttA 1 } + nA l {j, ttA 1 }) + 0(h n+2 ) 
= N l + h n+1 {j, F} + 0(h n+2 ) 

where we have used that (itA^^A 1 ) = P. 

We want to choose 7 so that {N n , H 3 }+ = 0(h n+3 ). Computing, we have 

(6.12) 

{N'\ iP'}* = {N l + h n+l { 1: P} , s(F 3 oN) + h n+1 ({G, h j } + g j o /))}* + 0(K n+3 ) 
= ift n+2 {{7, P}, h 3 } + ifr n+2 {P,{G : h 3 } + g 3 o 1} + 0(h n+3 ) 
= ih n+2 ({{7, P}, h j } + ih n+2 {P, {G, h j }) + 0(H n+3 ). 

where we have used the decomposition of theorem |2 Choosing 7 = G and using 
the Jacobi identity, we finish theorem [3J 

Remark. 7 is determined only up to a function of the actions, and thus the A 1 are 
not determined uniquely. To what does this freedom correspond on the quantum 
level? Consider the following transformation of the A 1 : 

(6.13) A 1 e -iGoN/H^i e iGoN/n^ 
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where G £ C°°(]R M )[[S]]. This transformation leaves the N* and the A algebra in- 
variant. On the level of wavefunctions, this transformation of the creation operators 
amounts to an N l dependent rephasing. 

Using the commutation relations, the transformation can also be written 

(6.14) A i -» A i e iGi0 t l l h , 

with Gi = G o (N + h/2) - G o (N + h/2 - he.,). At lowest order in h, we have 
s(Gi)/h = (diG) o I = {6\GoI}, and thus A i -> AV^< Go7 >. This is just the 
classical freedom in the origin of the angle variables. The indeterminacy of 7 is the 
manifestation in Eq. H6.13[l at higher order in %. 

Of course, Eq. . 1 1|1 shows the N l symbols themselves are uniquely determined 
at each order. Later, when we work with the Moyal star product, we will see directly 
that the 0(h 2 ) part of N l is uniquely determined from the H l . 

7. Example: Quantization of Moyal Symbols 

In this section, we will specialize to the Moyal star product and show how our 
prior analysis leads to higher order Bohr-Sommerfeld quantization rules for some 
quantum integrable systems of Weyl symbols^. 

The assumptions are now as follows: we have a set of M symbols H l which 
Moyal commute and whose principal symbols h 1 = nH 1 give rise to a global set 
of action angle variables {1,0)- We suppose that the series for H l are even in K. 
Further, we suppose that the functions z 1 = \j21 l cos 9 l and z t+M = — \J 2P sin 9 l 
are C°° and are Darboux coordinates. 

To achieve a quantization rule beyond the usual WKB order, we need to 

• extend the z 1 to symbols Z % satisfying irZ 1 = z l and {Z 1 ^^}^ — iKJ^ + 
0{h 5 ). 

• modify the Z l so that the associated number operators N l satisfy {ft 1 , iV- 7 }* = 
0{h 5 ). 

This is the work of section 17.21 and 17.31 respectively. The first problem will 
require some techniques. Before developing them, we must review the Moyal star 
product. 

7.1. The Moyal star product. The Moyal star product is a star product on 
T*(R M ) phase space, with the standard Poisson bracket 

(7.1) {f,9} = J ij d i fd j g = d i fd i g. 8 

Here, the partial derivatives di are with respect to the standard (x,p) coordinates 
on T*(R M ), and J is as before. The star product, for f,g£ C°°(]R 2M )[[fi]], is given 
by the Moyal formula 

(7-2) /*<? = J2m(t) 

where 

{f,g}o = f-9 

(7.3) {f,gh = {f,g} = d l fd l g 

{f,gh = d n d l Jd n d^g 



8 We use J to raise and lower indices. For example T l = J 1J Tj and Ti = JijT J . 



S'DARBOUX COORDINATES AND WKB APPROXIMATIONS IN DEFORMATION QUANTIZATION 



and so on. The Moyal commutator {/, g}± = f * g — g * / is an odd series in h and 
because of this, as we shall see, we can develop expansions for Z l , N 1 , etc as even 
series. Accordingly, the results from Sections [3] and will differ somewhat. 

The notation of Eq. (|7.3|l is not always convenient. To avoid proliferating indices, 
we will resume a diagrammatic notation introduced earlier in We define / — > 
g = {/,<?}, and, for more complicated diagrams, we exhibit an example conversion 
of an arrow to more explicit notation: 

D 

(7.4) = diDA - d l B »C 

A B C 



7.2. Construction of Moyal s'Darboux Coordinates. We need symbols Z l = 
z l + H 2 Z l 2 + . . . satisfying 

(7.5) {Z i ,Z j } ir = ihJ ij + 0(H 5 ). 
We have found such Z\. First define 9 

(7.6) r abc = / z b z c . 

Then 

(7.7) ^ = _Lr Qbc z^r abc 

satisfies Eq. lf?3)l . 10 

The key to verifying this will be the following identity, which we present without 
motivation: 



(7.8) {z d ,T abc } = 



9 The fact that the T abc are the connection coefficients for the flat symplectic connection in 
2* coordinates is not central to the derivation, but it is highly suggestive, especially in view of 
Fedosov's employ of a symplectic connection in constructing his star product. For now, we 
simply note that r ai>c is completely symmetric in its three indices. 

^Despite the known non-uniqueness of the Z2, one wonders if this construction is canonical 
in some sense. 
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The proof is as follows. Because the bracket {z , z } is constant, we have 



= 



+ 



= z d 



A + 



With rearrangement, the identity follows. 
A resolution of the arrow, 

(7.9) -^ = ^z a z a ^, 
will also prove useful. At order h 3 in {Z t , Z^}* we have 

(7.10) ih 3 Z{} - i <- i) - ±{z*, *} 3 
According to Eq. I|7.5|l . this must vanish. Working with the Z\ terms, we have 

{z\ Z{}-i~ j = ±{z\r ahc {zi,T abc }} - i <- j 

= 1 {{z\ T abc }{z\ r Qbc } + T abc {z\ {z\ T abc }}) - t <- j 

= ±{ z \r abc }{ z i,r abc }-i~j 



(7.11) 



-{z\T abc }{z\T^} 



1 

24 



We have used the Jacobi identity, the identity (|7.8(l . and the resolution of the 
arrow. Equation (|7.5|l is verified. 

We now have an explicit form for the first two terms in Z % = z l + h 2 Z\ + . . .. 
According to Theorem ^ the remaining terms in the Z % series can be completed so 
that the Z % are s'Darboux coordinates. To make the forthcoming analysis clearer, 
we shall now assume this has been done, although none of the results depend on 
the specific form of Z\, etc. 
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7.3. Good number operators. Associated to our Z l are number operators N % . 
They are easily expressed in terms of A 1 = a 1 + h 2 A\ + 0(h 4 ). 



N l — A 1 * A 1 

ih r _; h 2 

~2 



(7.12) = « V + « 4 } - y K, a'h + + a'Al) + 0(H 3 ) 



= p - \ - f {a*, a>} 2 ~ §r afcc r - + 0(h 3 ) 

The last two terms define 

(7.13) N\ = - l -{a\ a l } 2 ~ ^T abc P -> T abc . 

The appearance of an O(h) term 11 means that will be more convenient to develop 
F i as iT = F ; 5 (N l + h/2). 

Already, we have {iV J ,iP}* = 0(h 3 ). We would like to modify Z l so that 
{N\W}„ = 0(h 5 ). Taking F' 1 = f\ 12 we have W - F i o (N + h/2) = 0(h 2 ). 
By the second part of Theorem H H l — f l o (N + h/2) + h 2 G l for G l in which 
ixG % = {G, h 1 } + g l o I. The {G, h % } term in ttG 1 may be isolated subtracting the 

angle independent terms from H l — p o (N + h/2). Defining a complement to the 
angle average by )/( = /- (/), we have {G, h 1 } = W - f o (N + h/2)L 

Using a more accurate expansion of o, [^] we have 



(7.14) 

h 2 

16 

24 



fo(N + h/2) = h i - —{P,I k } 2 {d j d k f) o I 



(p -^i k ^- I'Wjdkdif) o / + tfN'idjf) o / + o{h A ) 



h l + h 2 K\ + h 2 n zj Ni, + o(h 4 ), 



where, for convenience, we have defined the frequency matrix f2jj = (Bjf 1 ) o / and 
(7.15) Kl = ^-{P j^ 2 { djdk p) o J + ^P -^P^I 1 (fydkdtf) o L 



We now have 

BG \ 

(7.16) {G, h*} = q^SI* = )Hi -K\- N'iUj \ 

and may solve for dG/dd k : 

BG 

(7-17) w = )V7 k \H k -K k ){-)Ni{. 



^This term is, of course, related to the Maslov index. 

12 Recall that /' is h? expressed as a function of its action variables, i.e. f l oI = h l 
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As in TheoremEl we define Z a = Z l + h 2 {-f, z l } + 0(h 4 ), 13 and choose 7 = G, so 
that {A n -kA n , Hi}* = 0(h 5 ). The new number operators are given by Eq. (|6.11(l : 

N !i + h/2 = N l + h/2 + h 2 {G, P} + 0(h 4 ) 

= r + h 2 N i 2 + h 2 ^ 

(7.18) 2 dP 

= P + h 2 N* + h 2 >Or fc 1 (^| - K k 2 )( - h 2 )N*{ 

= P + h 2 (Ni) + h 2 )nr k \H%-K*)( 

Note that (Ng) = (N l 2 ). 

Finally, we can obtain the quantization condition. We require 

H l = h l + h 2 H l 2 + 0(h 4 ) = F i o (N' + h/2) + o(h 4 ) 

(7.19) = f o (N' + h/2) + h 2 F* o (N' + h/2) + 0(h 4 ) 

= h i + h 2 K\ + h 2 a i:i N 2 3 + h 2 + h 2 F l 2 oi + o(h 4 ) 

By the work of Theorem |2 we know this equation is satisfied by taking 

Fi o I = (HI - Kl - n^N' 3 ) 

(7.20) 

= {Hi-K 2 --a ij N 2 1 ). 

F = /' + h 2 F 2 is the second order Bohr-Sommerfeld quantization rule for quan- 
tum integrable systems of Moyal symbols. Equation l|7.2L)[l contains two important 
terms: the K 2 comes from knowing the symbol of a function of an operator, while 
the N 2 term comes directly from the s'Darboux construction. 

8. Construction of s'Darboux coordinates in the Fedosov star 

PRODUCT 

Let M. be a manifold with symplectic form u> and a torsionless symplectic connec- 
tion V, i.e. one with Vw = 0. From {lo, V), Fedosov constructs a star product. |2] 
To summarize his result, given any / S C°°{M.) and Darboux local coordinates 

(n) 

z l , there are functions t which are completely symmetric in the numbered 

indices 14 such that Fedosov's star product is 

(8.1) 

f*9 = / (0 V 0) + § f^aS 1 + ■ ■ ■ + - (f ) " tfLy di • • • A + 

where w y is the Poisson tensor with respect to the z % . For early f^ n \ the results 
are 

/ (0) = / 
f (1) - V- f 
{ ' =V V- f 

Jii,i 2 ~ v 'i V '2J 



/i 1 ,i 2 ,i 3 ^ilVj 2 Vj 3 / Ri 1 i 2 i 3 3 V j f 

13 



Here we have one extra power of h 
4 For this section, we establish a com 
indices is a symmetrization over those indices 



14 For this section, we establish a convention that implicit any expression containing numbered 
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where the derivatives are taken with respect to the x % coordinates, and where 
Ru i 2 i 3 i — u)i x i R l i 2 i 3 k ^ k is the Riemann curvature tensor in those coordinates. 15 
Given Darboux coordinates z a , we wish to construct s'Darboux symbols Z a = 
z a + K 2 Z2- It is natural to guess that Z^ is as in Eq. (|7.7fl . but with T the new, 
possibly curved, symplectic connection written in the z a coordinates. We shall 
verify this. 

First, we need to prove a generalization of Eq. (|7.8(l : 

(8-3) V ai V a2 V a3 Z Raia2d3 — $ E aia2(l3 

where T ahc are the coefficients of the symplectic connection in z a coordinates. 

To prove this, we need only the fact that when r abc for a symplectic connection 
is written in Darboux coordinates, it is completely symmetric in its indices. [2] Let 
d a = d/dz a . First, we have 

V Ql V a2 V Q3 z d = V ai V a2 ^. 



— — V V d 

( 8 - 4 ) ~ ; 



0,3 
3«2 

8 r d + r CT v d + r CT r 

^Ctl 0302 1 aiQ3 (T«2 fli02 



0,30,2 0,10,3 (T02 ai02 &O3 

= —8 V d 4- 2Y a Y d 

0302 ' aia3 x aa2 ' 



where, in the last step, we have used the implicit symmetrization on numbered 
indices and the symmetry of T abc . Using the same property of r ahc , we also have 

p d a pd 1 ndp , pr; pd , p>) pd 

, , - rl -aia 2 Q3 — u a3 L ata-2 ' u L a l a 203 ~r L ai a 3 L T)a 2 ' a 2 a 3 L rjai 

— -3 v d 4- B d r 4- 2r ,? r d 

"03^ aia 2 T " 1 CHCI2O3 1 ■ 4J - aia 3 A 7702 ' 

Subtracting these two results, we obtain Eq. (|8.3|l . By a calculation very similar 
to the one in section 17^1 we can show that 

(8.6) Z d = z d + ^T abc {z d ,T abc } 

satisfy {Z\ Z*} = ihJ ij + 0(h 5 ). 

Following our previous development, we obtain the EBK formula in the Fedosov 
quantization by replacing the T in Eq. i|7.2U|l with this new V, and by changing the 
ordinary derivatives to covariant derivatives. 



9. Conclusion 

In this paper, we have developed a method for calculating spectra of functions 
which comprise a quantum integrable system under a star product. 

The methods developed here may have other applications. A future, more official, 
presentation may also include thoughts about 

• star product quantization rules on cohomologically non-trivial manifolds; 

• star Lie algebras from Poisson Lie algebras and applications to the semi- 
classics of spin; 

• Heisenberg evolution of s'Darboux coordinates. 



It seems Fedosov's own expression for /< 3 ) contains an error. He has 1/4 as the coefficient 
of R; we think it is just 1. 
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